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The critical indices a' , f3, 7' and 8 of the Quark Gluon Bags with Surface Tension Model with the 
tricritical endpoint are calculated as functions of the usual parameters of this model and two newly 
introduces parameters (indices). They are compared with the critical exponents of other models. 
It is shown that for the newly introduced indices x = and £ T < 1 there is a branch of solutions 
for which the critical exponents of the present model and the statistical multifragmentation model 
coincide, otherwise these models belong to different universality classes. It is shown that for realistic 
values of the parameter x the critical exponents a' , /3, 7' and 8 of simple liquids and 3-dimensional 
Ising model can be only described by the branch of solutions in which all indices except for a' 
correspond to their values within the statistical multifragmentation model. The scaling relations 
for the found critical exponents are verified and it is demonstrated that for the standard definition 
of the index a' the Fisher and Griffiths scaling inequalities are not fulfilled for some values of the 
model parameters, whereas the Liberman scaling inequality is always obeyed. Although it is shown 
that the specially defined index a' g recovers the scaling relations, another possibility, an existence 
of the non-Fisher universality classes, is also discussed. 
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I. INTRODUCTION 

Investigation of the properties of strongly interacting matter equation of state has become a focal point of modern 
nuclear physics of high energies. The low energy scan programs performed nowadays at CERN SPS and BNL RHIC are 
aimed at the discovery of the (tri)critical endpoint of the quantum chromodynamics (QCD) phase diagram. Despite 
many theoretical efforts neither an exact location nor the properties of the QCD (tri)critical endpoint are well known 
PP. Therefore, the thorough theoretical investigation of the QCD endpoint properties are required in order to clarify 
whether this endpoint is critical or tricritical. 

Here I would like to study the critical exponents of the QCD matter tricritical endpoint (triCEP). Since such a 
task, unfortunately, cannot be solved within the QCD itself and even the possibilities of the lattice QCD are nowadays 
very limited in this respect, then, unavoidably, one has to use some models of the triCEP. The most popular models 
of this kind are the quark-meson model (3j H] and the extended Nambu-Jona-Lazinio model [5] ■ Their triCEPs are 
generated by the intersection of the deconfinement phase transition (PT) line and the chiral symmetry restoration 
transition line. Since both kind of these models are the mean-field ones there is no reason to expect that their critical 
exponents could differ from that ones of the Van der Waals model equation of state [B]. Although very recently 
there appeared an interesting attempt to go beyond the mean- field approximation in the quark-meson model [7], 
the resulting temperature of the triCEP is about 32 MeV only, which is too low and unrealistic compared to the 
lattice QCD predictions [2 . Therefore, to study the properties of non-classical triCEPs one has to investigate the 
non-mean-field models. 

Recently there appeared a comprehensive analysis [S] of the QCD phase diagram with the triCEP within an extended 
version of the gas of bags model [5], which, however, is based on a physically inadequate assumption. Thus, in [5] 
it is assumed that the Fisher topological exponent t (denoted as a in 8J) of the Hagedorn mass spectrum should 
depend on the baryonic chemical potential. However, there is no reason to believe that this exponent characterizing 
the fractal properties of the surface of large bags should depend on chemical potential. Moreover, the results on the 
size distribution of clusters obtained for the 2- and 3-dimensional Ising model by the Complement method |10j clearly 
show that even far away from the critical point the Fisher exponent r has the same value as at the endpoint. 

Therefore, instead of model [8] here I consider another extension of the gas of bags model i,9, which is known as the 
Quark Gluon Bags with Surface Tension Model (QGBSTM) [TTJ [12] . This exactly solvable model employes the same 
mechanism of the triCEP generation which is typical for the liquid-gas PT and which is used in the Fisher droplet 
model (FDM) [T3] and in the statistical multifragmentation model (SMM) [T21 [TS] : the endpoint of the 1-st order PT 
appears due to vanishing of the surface tension coefficient at this point which leads to the indistinguishability between 
the liquid and gas phases. However, the surface tension coefficient in the QGBSTM has the region of negative values, 
which is a principally different feature of this model compared to the FDM, SMM and all other statistical models 
of the liquid-gas PT. Note that just this feature provides an existence of the cross-over at small values of baryonic 
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chemical potential li in the QGBSTM with triCEP [TT] and with CEP [T5], and also it generates an additional PT 
in the model with triCEP at large values of ix. Therefore, it is very important and interesting to study the critical 
indices of such a novel statistical model as the QGBSTM, to determine its class of universality and to examine how 
the latter is related to that ones of the FDM and SMM. 

The work is organized as follows. A brief description of the QGBSTM with the triCEP is given in Section In] In 



Section III the QGBSTM is analyzed in details and its critical exponents are calculated. Section [TV] is devoted to the 



analysis of scaling relations between the found critical exponents. Conclusions are given in Section [V] 



II. QUARK GLUON BAGS WITH SURFACE TENSION MODEL 

An exact solution of the QGBSTM was found in [11] . The relevant degrees of freedom in this model are the quark 
gluon plasma (QGP) bags and hadrons. The attraction between them is accounted like in the original statistical 
bootstrap model |17j via many sorts of the constituents, while the repulsion between them is introduced a la Van 
der Waals equation of state 0IUJ. An essential element of the QGBSTM is the T and ll dependence of its surface 
tension coefficient T£(Ts,yu) (here £(Ts,/i) is the reduced surface tension coefficient). Let us denote the nil line of 
the reduced surface tension coefficient as T^(li), i.e. £(Xs,/x) = 0. Note that for a given /x the surface tension is 
negative (positive) for T above (below) Is(/i) line. 

In the grand canonical ensemble the pressure of QGP and hadronic phase are, respectively, given by 

PQ{T,fi) = Ts q (T,li), (1) 
Ph(T,h) = T[F ff ( S H,T )M )+u(T, M )J T (A S ,£)], (2) 



where the following notations 



F h (s h ,T,li) = ^.g je ^-^0(T,m,), (3) 



oo 



7 T (A S ,E) = J ^e" A — (4) 
are used. Here sh = PH ^'^ , As = sh (T, fi) — sq (T, ll) and the particle density of a hadron of mass mj , baryonic 
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(p 2 J. m 2 -! 1 / 2 



charge bj, eigenvolume Vj and degeneracy <?j is denoted as <j)j{T,mj) = 7^2 Jp 2 dpe t . The QGBSTM is 



solved for a wide range of the functions u(T,li) and sq(T, //) [TT]. These functions are the parameters of the present 
model and here, as in [TT], it is assumed that the functions u(T, fx) and sq(T, li) and their first and second derivatives 
with respect to T and ll are finite everywhere at the T — ll plane. As it was shown in [11) . the 1-st order deconfinement 
PT exists for the Fisher exponent 1 < r < 2. In the continuous part of the spectrum of bags defined by Q the surface 
of a QGP bag of volume v is parameterized by the term v* . Usually one chooses x = | in 3-dimensional case (or 
x = for the dimension d) , but in what follows it is regarded as free parameter of the range < x < 1 . 

The system pressure, that corresponds to a dominant phase, is given by the largest value between pq and pu for 
each set of T and ll. As usual, the deconfinement PT occurs when pressure of QGP gets equal to that one of the 
hadron gas. Suppose that the necessary conditions for the deconfinement PT outlined in [TT] are satisfied and its 
transition temperature is given by the function T c {ll) for /i > fi cep . The necessary condition of the triCEP occurrence 
is that at the nil line of the surface tension coefficient there exists the surface tension induced PT of 2-nd (or higher 
order) [TT] for fx > [i cep and T c (ll) < T^(li) for these li values. The both of these inequalities become equalities only 
at triCEP. Moreover, at triCEP the phase coexistence curve T c (^) is a tangent (not intersecting!) line to the nil line 
of the surface tension coefficient Ts (/i) . For li < LL cep the deconfinement PT degenerates into a cross-over since in this 
region £(/i) < [TT] and, hence, the system pressure is defined by a solution of Eq. A schematic picture of the 
phase diagram in li — T plane is shown in Fig. [T] 

An actual parameterization of the reduced surface tension coefficient S(T, ll) is taken from |11| (£ = const): 



sign(T E ( M )-T) , (5) 



but here it is written in a more general form which allows one to perform the calculations for the integer and real 
values of power In what follows the coefficient <tq is assumed to be a positive constant, i.e. ctq = const > 0, but 
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FIG. 1: [Color online] A schematic phase diagram in the plane of baryonic chemical potential /i and temperature T. The 
dashed curve indicates the nil line of the surface tension coefficient Ts(/x), below (above) which the surface tension coefficient 
is positive (negative). The deconfinement PT line T c (/j,) is shown by the full curve for fi > /i cep whereas the surface tension 
induced PT is depicted by the long dashed curve for [i > /i cep . A cross-over (shown by the short dashed curve) takes place 
along the line Te(^i) for \i < fi C e P - The cross-over and PT regions are separated by triCEP (filled circle). At this point the 
curves T c (fi) and TsQi) are tangent to each other. 



it is easy to show that the obtained results hold, if <jq > is a smooth function of T and /i. Here it is appropriate 
to say a few words about the negative values of E. Note that the existence of regions with E < is a distinctive 
feature of QGBSTM compared to other models. There is nothing wrong or unphysical with the negative values of 
surface tension coefficient, since TE))" is the surface free energy of the bag of mean volume v and, hence, as any free 
energy, it contains the energy part e sur f and the entropy part s sur f multiplied by temperature T [13] . Therefore, at 
low temperatures the energy part dominates and surface free energy is positive, whereas at high temperatures the 
number of bag configurations with large surface drastically increases and it exceeds the Boltzmann suppression and, 
hence, the surface free energy becomes negative since s sur f > e "^f . Such a behavior of the surface free energy can 
be derived within the exactly solvable model of surface deformations known as Hills and Dales Model [TH] . Moreover, 
very recently using the relation between the tension of confining color string and the surface tension of the QGP bag 
derived in [TH] it was possible to demonstrate [HJ [20] that at the cross-over region the surface tension coefficient of 
large bags is unavoidably negative. In addition this approach allows one to determine such an important parameter of 
the QGBSTM as the value of triCEP/CEP temperature T cep = 153.9 ±4.5 MeV [5T] under the plausible assumptions 
which are typical for the liquid-gas PT [13l [18] . Note that this value of the triCEP / CEP temperature is in a good 
agreement with the lattice QCD results discussed in [3J. 

In the vicinity of triCEP the behavior of both the deconfinement PT curve and the nil surface tension coefficient 
line in the /x — T plane is parameterized via a single parameter £ T > 0: 



T, 



cep 



Ts(fi) ~ T cep - T c (p) ~ (fj, - /i cep ) 4 



(6) 



since, as discussed above, T c ([i) and Ts(//) are tangent to each other at triCEP. This is a new index which was not 
considered both in the FDM and SMM. As it will be shown below this index is responsible for a new universality class 
compared to other exactly solvable models. 



At the phase coexistence curve the Clapeyron-Clausius equation can be written as 



The 



T=T C 



Sh — Sq 
Ph-Pq 

entropy density Sa aud the baryonic density pa are, respectively, defined as T and /i partial derivatives of the 
corresponding pressure pA, where A 6 {H, Q}. Then using and ^ the Clapeyron-Clausius equation can be 
explicitly rewritten as 



dfi c 



9Fh 
dT 



OFh 



£/t(0,£) 



dT 



( 9F H 
\ ds 



1 



+ 



^J T (0,S) + ^ 



dfj, J 



dfj, V ds 



as 

dT 



§ U I T _„(0,E) 



(7) 



T=T C 
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Let us parameterize the behavior of the numerator and denominator in ^ at the triCEP vicinity as 

-l) ~ (T cep -T c (p)Y 

- I | ~ {T cep - T e {n))x, 

T=T C 



\dF H 


H^ir(0,E) + 


5s Q 


(dF H 


dT 


<9T 


^ ds 


dFjj 


h ^/ T (0, £) + 


3sq 


(dF H 


dfi 


dp 


\ ds 



(8) 
(9) 



where x > max(0, 1 — ) denotes another new index. The latter inequality follows from the fact that the integral I T 
and functions Fh, u, sq together with their derivatives are finite for any finite values of T and p. It is necessary to 
emphasize that the index \ unavoidably appears from an inspection of the Clapeyron-Clausius equation. Since the 
latter is a direct consequence of the Gibbs criterion for phase equilibrium, then an introduction of the parameter x is 
quite general and could be done for any model with the PT of the liquid-gas type. So far, the index \ was never used 
for the calculation of critical exponents. For example, in the SMM case this parameter is set to zero by construction 
[221 [231. 



III. THE CRITICAL INDICES OF THE QGBSTM 



The standard set of critical exponents a', f3 and 7 [51 [23] describes the T-dependence of the system near triCEP: 



Ap 
AK T 



\t\~ a , for t < and p 
\tf, for t < 0, 
|t|- 7 ', for t < 0, 



Pcepi 



(10) 

(ii) 

(12) 



where Ap = (pg — Ph)t=t c defines the order parameter, C p = ^(ffOp denotes the specific heat at the critical density 
and AKt = (K^ — K®)t=t c is the discontinuity in the isothermal compressibility Kt = ^(§^)t across the PT line, 

. The critical isotherm shape is given by the index 8 [51 12"4"] 

(p cep - for t = 0. (13) 



the variable t is the reduced temperature t 



T-T c , 

Veep P r 



Hereafter the tilde indicates that T = T cep . 

The calculation of a' requires the knowledge of the specific heat behavior C p = ^(f^) P along the critical isochore 
P = Peep inside the mixed quark-gluon-hadron phase. In contrast to the works [221 [24], the famous Yang- Yang formula 
P5] is not used here to calculate C p , since this formula leads to a less convenient representation. Therefore, below 
there are given some details of the specific heat evaluation. 

As usual the entropy density and baryonic density of the mixed phase are defined via that ones of the pure phases 
and the parameter A, which is the volume fraction of hadronic phase (i.e. the volume fraction of QGP is, respectively, 
1-A): 



P\T=T C 
S\t=T c 



\Ph\t=t c 
\Sh\t=t 



(1 - A)pq| t =t c , 
(1 - \)Sq\t=T c 



(14) 
(15) 



Varying A from to 1 one can describe any state in the mixed phase for fixed p or T, since these variables are not 
independent inside the mixed phase. Replacing p\t=t c by p cep in Eq . (|1 4[), one can calculate the total T-derivative 
of A along the critical isochore p = p cep - Then for p = p cep from Eq. (|l5[) one finds 



T 

Pcep 



A 



dT 



(Sh — Sq)t=t 



dp c d 

It It 



(ph - Pq)t=t c 



( d(S Q \ T =T c ) dp c d{p Q \ T =T c ) 
\ dT dT dT 



Using the Clapeyron-Clausius equation it is possible to rewrite the specific heat ( 16 1 in the form 

d 2 p c 



T 

Pcep 



(PQ - Peep)T=T c 



dT 2 



( d(S Q \ T =T c ) dp c d{p Q \ T=Ta ) 



dT 



dT 



dT 



(16) 



(17) 



where the definition ( 14 1 for A was used. 



The critical exponent a' describes the temperature behavior of the most singular term in Eq. (17 1. Expanding 



Pq\t=t into series of i-powers and using the parametrization mh of the coexistence curve one can show that the first 
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term in |l7| behaves as f 1 ™ 1 - 1, s T ' + 5 T 2 . Since the entropy density and baryonic density of QGP are defined via the T 
and fj, derivatives of the function Tsq(T, /z), which is a regular function together with its first and second derivatives, 

then the singularity in the second term of (17) appears from (^fr^J ~ t( T 2 . Clearly, the second term in (17) is 

a non- vanishing constant, if T-derivative of /J. C (T) is finite at triCEP. Accounting for this fact, one gets the critical 
exponent a' as 



a' = 2 - 2 min 1, 



e 



(18) 



This equation shows that a' > for £ T > 1 only, otherwise a' — 0. As it will be seen below only the index a' > 
depends on £ T , whereas other critical exponents always have the branch of solutions which is independent on £ T . 
Compared to the FDM and SMM, such a property leads to a principally new and unique possibility to choose the 
index a' independently of other critical indices. 

To calculate the critical exponents (3 and 7' it is necessary to analyze the behavior of the integral 7 T _ g (0, E) for 
small positive values of £, which is the case, when T — T c (fi) approaches T cep . In the limit E — > +0 the integral 



I T - q (0 7 E) remains finite for r > 1 
substitution z e wSi yields [2"2l l2"i 



q and diverges otherwise. For r = 1 + q it diverges logarithmically. Indeed, the 



/ T -,(0,E) = E 2 



dz 



z T-q 



(19) 



The condition r < q + 1 guaranties a convergence of the integral (19) at its lower limit. The finite values of 



J T _q(0, E) at the upper limit is provided by an exponential factor. Summarizing these results for E — > +0 one can 
write 



Jr-,(0,E) 



E T = " , for r < 1 



InE, 
const, 



for r = 1 + q 
for r > 1 + q 



E min(0,^) 



(20) 



where the limit {x v lncc) x _>.o ~ % v w &s used. 

To calculate the index 5 one should investigate the integral Q in the limit As — > +0 and E — > —0, but for a 
— > —0 which, as shown in [11], is the case (also see Eq. (28) for details). Then one obtains 



condition 

As 



7 T _ g (As,E) 




As i: 



i(0,T-l-g) 



(21) 



Using the definition of the baryonic density and Eqs. ([T]) and ([2| one can find the baryonic density discontinuity 
across the dcconfincmcnt PT line: 



Ap = T c 



dF H 



§M0,£) 



dF„ 
ds 



1) 



Ulr-^O, E) 



ds 



•«I T _ 1 (0,S) 



T=T C 



According to the parameterizations ( 5 ) and |6| one obtains that along the PT line E ~ and |p 
hence, the temperature dependence o: 

r) + \ 



(22) 



and, 



Ap in ( p2| ) can be straightforwardly found from Eqs. (|9| and (20). This yields 



x 



£ min(x, r - 1) - |V 



for 
for 



X < ^min(x,r- 1) - ^ 
X> i min(x, r - 1) - i 



(23) 



The structure of this equation is obvious: its first term corresponds to the denominator of Eq. (22), whereas the 



second one describes the leading term of the numerator of ( 22 ) 



To find the critical exponent 7' one has to calculate the isothermal compressibility Kt = ^(§^)t both for QGP and 



EE if^ 
p ^ dp 

for hadronic phase. With the help of the baryonic density definition one can rewrite the isothermal compressibility 

1 p& 

as Kt — \ n-Tf ■ This result clearly demonstrates that the QGP contribution into the AKt is negligibly small since 



the QGP pressure is given by the function sq(T,h), which, according to the QGMSTM assumptions [TT], is finite 



G 



together with its first and second derivatives for all finite values of T and p. Therefore, at triCEP AKt can diverge 
due to the hadronic phase contribution only, i.e. near triCEP AKt — K^ ■ Therefore, using Eq. ^ and keeping at 
triCEP the most singular terms one finds 



AKj 



^l-^+u[ M (0,E) 



J r - 2 (0,S) _ fdAs" 
J T _i(0,S) V. <)// / T - V ; ( ,,, 



C 



y = - - 2/3. 



(24) 



T=T C ( M ) 

From the definition of As it follows that along the PT line = - t p . Then from Eqs. (po|, |5|, (fsj) one has 



(25) 



At the critical isotherm one can use the definition of As to get 



P-Pcep =T cep (As + . 



Q ~ s Q\cep) 



T cep As + Afi 



dp 



(26) 



where in the second step one has to expand sq in powers of Afi = ji — p cep and keep the linear term. Similarly one 
determines the deviation of the baryonic density taken at the critical isotherm that is lying outside the mixed phase 
from the baryonic density at triCEP: 



P Peep T c , 



OAs Osq_ _ ds Q 
dp dp dp 



cep i 



dAs A d 2 s c 
~dp + Afi — 



dp 2 



(27) 



cep j 



Now it is clear that the behavior of As near triCEP should be analyzed in order to calculate the critical exponent 
S. Let us consider the case \ = first. Such a case is typical for the SMM [TSJ [55] . Since at triCEP As = 0, then 

substituting the expansion I T {As, E) « J T (0, 0) - As/ T _i(^, §) - EJ T _ X (^, f) [TT] into Eq. Q for the cross-over 
states one gets 



As 



^(ff + f^(0,0) + ^(^ 



1 dp, ^ al r-l\ 2 > 2 > 



(28) 



where it is sufficient to keep only the first order terms of expansion, whereas for the case \ > one has to keep the 
second order terms as well. With the help of Eq. (21 1 one concludes that 



As 



A^ Jffl «(T-i,«) ) for 



z 1 C 



< 



<(t — — r — 1 



Ap^ 

which allows one to find the index S for \ = 



for 



> 



(29) 



max(r — 1,jht) — r — 1 



5 lx=o 



max(r-l,x) 

r-1 
2-t ' 



for 1 < f c , , < -Xr 



for 



max(r— l,>f) 



(30) 



> 



Note that in this case the inequalities | < r < 2 providing the existence of the 2-nd order PT at triCEP [TT] also 
guaranty the fulfillment of the condition S > 1. 

The result for the index S in the case x > can be found similarly: 



max(r — 1, k) 



1 



for 1 < 



max(r-l,jf] 



< 



(31) 



Since in some aspects the QGBSTM is similar to the FDM and SMM it is interesting to compare its critical 
exponents with that ones of the FDM [13] and SMM |22|. However, one has also to remember that in contrast to 
the FDM, both the QGBSTM and the SMM have a non-vanishing excluded volume of the constituents and the same 
range of the Fisher exponent 1 < r < 2 for the triCEP existence, whereas the FDM can be formulated for r > 2 
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only. Since the FDM and SMM implicitly treat the parameter x = 0, then it is most natural to compare their critical 



exponents with the QGBSTM results just for such a case. Eq. (30) shows that in the QGBSTM there is a regime, 
when its index 6\ x —o matches the SMM result 22] . Moreover, it is easy to see that in this regime all other indices the 
QGBSTM and SMM coincide for £ T < 1 

a' = 0, /?=i(2-r), r/ = ^( T -f) and 6\ x=0 = , (32) 

X X \ Z J Z — T 

although the present model has entirely new regime for £ T > 1, for which all its indices except for a 1 are equal to the 
corresponding exponents of the SMM, whereas a' can be chosen freely provided that the index £ T exceeds the value 
max(r-i^) ^ Q n ther hand the QGBSTM always belongs to a different universality class 1 < r < 2 (except for 
a special case r = 2), than that one of FDM in which r > 2. Thus, the spectrum of values of the QGBSTM critical 
indices is more rich than the corresponding spectra of the FDM and SMM since this model contains two new indices 
£ T and \- 

Let us demonstrate this using a few sets of critical indices. Since the present model contains five parameters, r, x, 
£ T and x, there exists an infinite number of possibilities to describe the four standard critical exponents, a' , f3, 7' 
and S. Thus, the critical exponents of the 2-dimensional Ising model [29] shown in the first row of the Table I can be 
exactly reproduced by many sets of parameters with the vanishing value of index \ ( see the nrs t row m the Table II) . 
It is interesting to mention that, if one uses the linear T-dependence of the surface tension coefficient ^ employed in 
the FDM for T < T cep , i.e. fixes £ = 1, then for the 2-dimensional Ising model one obtains x=\, which is typical for 
the dimension d = 2. However, in contrast to the FDM, where t = || > 2, this set of critical exponents is described 
by the value r = || < 2. 

It is interesting to analyze the critical indices of the simple liquids since near the deconfinement region the QGP 
behaves as a strongly interacting liquid pQ . Taking the wide range of values for the critical exponents of simple liquids 
[2"9"] (see the Table I) one can describe them in many ways. The set A in the Table II demonstrates one of such 
possibilities, but its value of the parameter x = (t — 1) < ~ 0.6069 is essentially smaller than the expected for 
3-dimensions value x sa |. At the first glance it seems that the accuracy of above 10% for the set A value of x which 
is obtained for r < 2 is acceptable, but a careful study of the surface free energy in various cluster models [121 [Ml I3"U] 
shows that for 3-dimcnsional case the value x = | holds with much better accuracy of about 4%. However, if one 
requires that the critical exponents of simple liquids are reproduced for x > 0.49, i.e. in such a way that both 2- and 
3-dimensional values of parameter k can be included, then it can be shown that the only existing solution corresponds 



to the SMM values for indices /3, 7' and <5 ( 32 ) with x — 0, while the index a' is fixed to its experimental value (see the 
set B for liquids in the Table II). In fact, the same outcome is obtained for the critical exponents of the 3-dimensional 
Ising model (see the Table II) which were found with very high accuracy in [31j and are given in the Table I. 

There are two important consequences that follow from these results. First, since the QGBSTM indices (3, 7' and 
6\x=o match that ones of the SMM and the latter do not depend on the parameter £ T , then there exist two relations 
for index r [22] 



' and t = 2- 13 . (33) 



1 + 5 j 1 + 2(3 

The first of these equalities follows from the definition of index <5| x= o = 537 > whereas the second one can be obtained 



directly from Eq. (32). Second, from (33) it follows that the index t ~ 1.826 ± 0.02 has a very narrow range for 
the set B of simple liquids and for the 3-dimensionl Ising model (see the Table II). Both of these consequences are 
important for QCD since its universality class is expected to match the class of the 3-dimensionl Ising model [2, 32, 33 . 
Therefore, on the basis of above results one can expect [H] that at triCEP the volume distribution of large QGP bags 
has a power like form v~ T with r ps 1.826 ± 0.02. The same conclusion for the mass distribution results from the fact 
that in the QGBSTM [11] the mass of large bags and its volume v > Vq are proportional to each other. 



IV. THE SCALING RELATIONS OF THE QGBSTM 

The well known exponent inequalities were proven for real gases by 

Fisher [26] : oi + 2$ + i > 2, (34) 
Griffiths [27] : a' + /3(l + 6)>2, (35) 
Libcrman [28] : 7' + /3(1 - 5) > 0. (36) 





a 


P 


7 


5 


2d Ising model 





i 
s 


7 

4 


15 


Simple liquids 


0.09-0.11 


0.32-0.35 


1.2-1.3 


4.2-4.8 


3d Ising model 


0.1096 ±0.0005 


0.3265 ± 0.0001 


1.2373 ±0.0002 


4.7893 ± 0.0008 



TABLE I: The critical indices of simple liquids [29] . 2-dimensional Ising model [29] and 3-dimensional Ising model 131] . 





X 


e 


T 


XT 


c 


2d Ising model 







31 
1C 


min(2x, f ) > ^ 


2x 


Simple liquids A 


0< X < 4 


20 
19 


20 , 23 
11 t X 44 


( T - 1)«7 

V -V720 


23^ 


Simple liquids B 





1.0526 ± 0.0055 


1.8255 ±0.0212 


x > 0.4947 


x- (1.92 ±0.026) 


3d Ising model 





1.0579 ±0.000055 


1.8272 ±0.000048 


x > 0.4999 


x- (1.8903 ± 0.00007) 



TABLE II: The QGBSTM parameters that describes the corresponding exponents given in the Table I. 



The corresponding exponent inequalities for magnetic systems are often called Rushbrooke, Griffiths and Widom 
inequalities, respectively. These inequalities are traditionally believed to play a fundamental role in the modern 
theory of critical phenomena. However, the real situation with the scaling inequalities (34)-(36) is not that trivial 
as it is often presented in the textbooks. Thus, as one can see from the Table III the inequalities (34)-(36) exactly 



hold for the Onsager solution only, whereas for simple liquids only the Fisher relation is obeyed and even for highly 
accurate numerical evaluation of the critical exponents of the 3-dimensional Ising model there are some problems 
with the inequalities (34) and (36). Moreover, a long time ago it was found [24 that the traditional definition of the 



exponent a' given by (10) may lead to somewhat smaller value than 2 staying on the right hand side of Eqs. (34) 
and ( 35 1 . A similar result was analytically found for the SMM 
SMM parameters [141 H~5] the right hand side of inequalities (|34 



3|j_which shows that for the standard set of the 
and (35) should be replaced by Therefore, it is 
interesting to verify the scaling inequalities for the QGBSTM indices ODlained here. 

Despite the usual expectations, the QGBSTM critical exponents do not obey the traditional scaling relations in 
general. Again, as in the SMM case, the Fisher and Griffiths inequalities are not always fulfilled, whereas the Liberman 
inequality i s fu lfilled for any values of the model parameters. Indeed, let's de monstrate the validity of the Liberman 
inequality (36) first. For simplicity, consider the case S\ x=0 = §Ef of Eq. (30), which is realized for x 
mentioned in the preceding section the QGBSTM indices f3, 7' and 6\ x —q for t 
exponents of the SMM and, hence, as in the SMM case [52], the Liberman inequality is fulfilled within the present 
model for any choice of ol . This, however, can be shown from the explicit expressions for the indices /3, 7' and <5| x — 0, 
i.e. from Eqs. (23), (25) and (30). Using these equations one obtains 



0. As it was 
lis case coincide with the corresponding 



lin ^0, min(r — 1, x) — ^ 



>0, 



(37) 



where the validity of the right hand side of (37) easily follows now from the inequalities min(0, . ..) < and r < 2. 
The Liberman relation analysis for other values of the index S gives the same result. Using the Liberman inequality 
and the explicit expressions for the QGBSTM critical exponents one can get the following result for the Fisher and 
Griffiths inequalities 



+ (3(6+1) <a' + 2/3 + 7' = 2 + 



- - nun i 1, 



(38) 



which holds for any value of the index x- This equation clearly demonstrates that the Fisher and Griffiths inequalities 
are not obeyed for the values of parameters satisfying the inequality J < 2min(l, Jj?). Moreover, a fulfillment of the 
Fisher scaling inequality does not guaranty that the Griffiths one is obeyed. 

In order to 'save' the scaling inequalities ( 34 ) and ( 35 ) it was suggested [53] to replace the index a' by a' s , where the 
index a' s describes the temperature dependence of the specific heat difference AC = (C PH — C pq )t=t c for two phases. 
Since at the coexistence curve the specific heat of each phase is defined by the total T derivative of the entropy, then, 
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Fisher: a' + 2/3 ± 7' 


Griffiths: a' + p(8 + 1) 


Liberman: 7' + /3(1 — 5) 


2D Ising model 


2 


2 





Simple liquids 


2.02 ± 0.0055 


1.9425 ±0.0055 


0.0775 ±0.0212 


3D Ising model 


1.99996 ± 0.00007 


2.000412 ± 0.005 


-0.000052 ± 0.002 



TABLE III: Scaling relations between the critical exponents taken from the Table [I] The uncertainties were calculated from 
their values given in the Table |I] using the error determination method for indirect measurements |34], 



using the Clapeyron-Clausius equation one can find AC as 



T c d 
Ph dT 



Pq ~ Ph dS Q 



PHPQ 



dT 



(39) 



where the entropy density Sq and the baryonic densities pq, pu are calculated along the deconfinement PT line. 
Then from the parameterization ^ and the definition of index /3 one gets 




1 



for <2 
for jr>2 



max 2 



1 

F 



(3 



1 



(40) 



Note that a' s > for Jr < 2 — j3 only. Using a' s from ( 40 ) instead of a' in the Fisher inequality, one can write 



a' s + 2(3 ± 7' = max ( 2, ^ 



c 



c 



,TJ \ H - P - £T ) ^ 2 - T. mm (°- 1 + ^ - T) > 2, 



(41) 



where the final result follows from the obvious inequalities max(2, ...) > 2 and min(0, ...) < 0. Similarly, one can prove 
the validity of the Griffiths inequality for the index a' s . Thus, the Fisher hypothesis [24] to replace a' by a' s recovers 
the scaling inequalities for the critical exponents, however, it does not seem that such a suggestion is the final solution 
of this problem. 



CONCLUSIONS 



The critical indices of the QGBSTM with the triCEP are determined and compared with the critical exponents 
of other models. The QGBSTM critical exponents are expressed in terms of the model parameters t, x, C an d two 
newly introduced indices £ T > and \ ^ max(0, 1 — The index £ T in |6| characterizes the behavior of the PT 
curve in the vicinity of triCEP in the plane of baryonic chemical potential and temperature, whereas the index x i n 
Q describes the temperature dependence of a certain combination of the /i-derivatives of the model spectrum in the 
same vicinity. 

Since in the FDM and SMM the index \ is implicitly set to zero, while the index £ T either docs not exists (FDM) 
or is fixed by other model parameters (SMM), then the spectrum of the values of critical exponents of the present 
model in more rich compared to those models. For the case x = and £ T < the QGBSTM reproduces the critical 
exponents of the SMM with triCEP, whereas for other choice of parameters x an d £ T these models belong to different 
classes of universality. Also the universality classes of the FDM and QGBSTM are different since the range of values of 
their index r is different (except for a singular case r = 2): r > 2 in the FDM and 1 < r < 2 in the QGBSTM. A very 
important result of the present work is that, if one requires that the parameter u describing the surface dependence 
on the bag volume should have the values typical for dimensions 2 and 3, then the critical indices of simple liquids and 
3-dimensional Ising model can be described only by the SMM expressions for j3, 7' and S exponents while the index 
a' differs from the SMM value. For these sets of critical exponents a very narrow range of parameter r = 1.826 ± 0.02 
is found. Such a prediction might be important for experimental searches of the QCD phase diagram endpoint since 
just this exponent describes the power law in the volume distribution of large bags at triCEP. 

The direct calculations show that for the standard definition of the critical index a' (found along the critical 
isochore) the Fisher and Griffiths scaling inequalities are not always fulfilled, whereas the Liberman inequality is 
obeyed for any values of the model parameters. In contrast to the SMM, in which the critical isochore belongs to 
the boundary of the mixed and liquid phases, the critical isochore of the present model is located inside the mixed 
phase and, hence, the conditions of the Fisher theorem [26] proving the validity of (34) are formally satisfied, but 
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the Fisher and Griffiths scaling inequalities are not fulfilled. Therefore, it is quite possible that instead of the Fisher 
suggestion one should search for an alternative solution of this problem and, thus, to admit an existence of other, 
non-Fisher, universality classes of critical exponents for which the right hand side of scaling inequalities should be 
modified. Hopefully, further theoretical and experimental studies of this problem will find its final solution. 

Acknowledgments. The author appreciates the stimulating and fruitful discussions with K. A. Bugaev. 
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